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THE MODEL OF THE COMPOSITIOU OF THE MAETIM ATMOSPHERE 


M. N.' Izakov, O/'F.'.TCrasitskiy 
Introduction 

A knowledge of tlie composition of the planet !s .atmosphere, .along 
with its thermal conditions and dynamics is necessary to understand 
the structure of the atmosphere and the processes that take place in it, 
as well as to solye the problems of the evolution of the atmosphere. 

During the past few years we have noted considerable progress in the 
study of the Martian atmosphere, including its composition, to a 
significant degree because of research using spacecraft. (See, for 
example, reviews [1-U] where there is a detailed bibliography). ‘ 

The Martian atmosphere (the average pressure at the planet's surface 
is 6 millihars, the average temperature is 220°K [1-5] ) contains approximately 
95/^ 2-3^ and 1=2^ Ar [ 6 - 7 ] . A large quantity of argon, obtained 

by indirect means [8]., is either excessive or belongs to some kinds-io'f 
local conditions. In addition, the atmosphere cbntains"rsnialihGompanents 
4the'-Trolumetriic 'Content telativettooGO^) : molecular oxygen 0^ (0.35“1.2)10-^ 

carbon monoxide CO- (0.ii-3.0)l0 "The content of water vapor in the 
atmosphere changes with the seasons and the latitude within the limits 
of approximately 5 to 50 p of the precipitated water, and sometimes O’-yj 
over the polar caps it falls to a fraction of a U, andiin lower lati- 
tudes it sometimes reaches 100 y [13-lT] • In the Martian atmosphere 
ozone is found; according to data of the tiltra-violet spectrometers on 
the Mariner 6, T and 9 spacecrafts in the low latitudes, the level of 0^ 
was less than the sensitivity level of the device, equal to 3 y-atm. 



(1 11 -atiu. — . 2. 687 • lO^^cra ^), and in the polar areas the level of O 3 
varied from 3 to 15 y-atm- , and sometimes reached 60 y-atm. [18—19] ; 
According to similar measurements on the "Mars-5" spacecraft, it was 
determined that in low latitudes the content of 0^ is not less than 3 y-atm. 
and that the maximum of the layer of 0 ^ is located at approximately 
40 kilometers [20]. The concentration of atomic hydrogen ate 
altitudes of 200-250 km (measured by the luminescence of the atmosphere 
in the lyman-a3.pha line of hydrogen 1216 Angstroms) is equal to 
n^(250) = (2-3) • 10*^cm~^ [21—24]. From the emissions of the Martian 
atmospl^ere in the lines of atomic oxygen we found that! at an altitude 
of 135 km (close to the ionospheric maximum) the concentration of 
the atomic oxygen in relationship to the concentration of CO 2 
is equal to n (135)-0. 5-352 n [25-28]. From the emissions in the 
CO zones we made the following evaluation: u^q (135) = 0. 3-1.0% n^Q 

2 

[26] . 

Theoretical models of the composition of the Martian atmosphere, 
permit in the first place, the coordination of different experimental 
data; secondly, prediction of the presence and altitude profiles of several 
atmospheric components that still have not been found in the experiments; 
thirdly, explanation of the surprising fact of the predominance of CO 2 
to high altitudes in the presence of its quick photodissociation and 
slow association of CO and 0 [29—32]. However, several hypotheses 
formulated when constructing these models need to be refined, not all 
of the computed concentrations concur with the experimental data 
and have significant differences between the data from different models. 

In this report we suggest a model of average global distributions 
of the atmospheric components from the surface tooan altitude of 250 km,. 
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in which we used a more accurate turbulent mixing coefficient, depending 
on the altitude; we have refined the rates of several chemical reactions 
and the rates of photodissociation of several components according /5. 

to the latest experimental data; we used data obtained by us earlier on 
the thermal condition of the upper atmosphere; we developed a detailed 
method of computations. ite a result we have a good agreement between 
the data of the models and all of the experimental data; we have explained 
possible variations in the concentrations of atmospheric components at 
different altitudes depending on the possible variations (long periods, 
for example - seasonal), the temperature, the humidity and the turbulent 
mixing; we have defined the relative role of the different processes 
determining the atmospheric composition of Marss Limitations with long 
period variations are determined by the fact that we have found 
steady-state distribution of the concentrations under global-average 
conditions. However, since these distributions are determined by 
integration in time, this subsequently permits the use of the model 
apparatus after modernizing it for studying both short-term (daily) 
and longer period (evolutionary variations in the Martian atmospheric 
composition. 

1. Construction of the Model 

Spatial distributions and temporary variations of atmospheric components 
are described in equations of continuity for the concentrations of the 
components which include the basic processes determining the composition 
of the process - the origination and destruction of particles in chemical 
and photochemical reactions and their transfer by molecular diffusion and 
turbulent mixing (see, for example, [4]). The composition is of course 
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connected with the thermal conditions and movementSi’in the atmosphere, 
so that a strict examination requires simultaneous solving of the 
continuity equations with the-;equations of movement and energy 
balance. However, in view of the fact that the basic component 
in the Martian atmosphere is CO^ which prei;ails .to high altitudes, 
we can initially find the thermal conditions and movements for the 
atmosphere from (as it was done in [33]) and, using the temperature /6 . 
profiles found there, we can solve the system of continuity equations for 
the components. 

An evaluation of the characteristic times of the different processes 
which* determine the planet’s atmosphere show that the distribution 
of the components (Z) at the altitude are basically determined hy 

K 

slower processes of transfer along the vertical, whereas the quick 
transfer along the horizontal (which is a result of the quasi-horizon- 
tality of atmospheric movements) to a significant degree decreases the daily 
variations. This fact allows us to ignore the derivatives along 
the horizontal in the transfer member and, thereby, to reduce the 
problem to a spatial— unidimensional one, having written the equation 

f 

of continuity for concentrations of the components'll^ in the following 
form: 
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Here, t - time, Z - tHe altitude; - the diffusion coefficient 


of the particles of the K-component, K, -’the coefficient of the 
turbulent mixing; T is the temperature; is the molecular weight 

of the k-component, y is the awerage molecular weight (il=n , n=En, ) 

X 

a is the thermodiffusion factor (significant only for hydrogen, 
ic 

'= 0.38); R is the universal gas constant; g is the acceleration 
' ° 

of the force of gravity; and are the speeds of the appearance 
and disappearance of the particles of the k-component in the 
photochemical and chemical reactions. 

The members describing the photodissociation of the molecules 
of the X.-type take on the form (see, for example, ’[:U, 3^1): 


I I'v ^ ^ ^ 


. I 


•), 


where A, V are the wave length and the photon frequency, X is the 

zenith angle, of the> sun, h is the Planck constant, .4.1, 4. n 

“ - ’ ’ A™ is the spectral 

stream of the solar ultraviolet radiation on the boundary of the 

atmosphere (the values of which are used in accordance with [35] , for 

X<1000 S, and according to [36], for large X); CJ^ is the absorption 

coefficient by the 51-molecules, is the photodissociation coefficient, 

equal to a^= when X 1000 S. and a^= when X < 1000 S 

(here b^ is the ionization coefficient); the values of b are 

X, bU2 

according to [37-39], of b^ according to [36, 40, 41], Q,. according 
to [36] , according to [42] , according to [43] , 



according to [44] . 

Tlie cross-section of the attenuation of radiation with 
carhon dioxide gas quietly decreases the wave length X 

up to X = 2050 %, and with larger X it decreases much more slowly 
[38] . iShemansky [38] showed that this deviation on the a^^^(X) 
curve is explained by the fact that here a cross-section of Rayleigh 
scattering is added to the cross-section of the absorption 'of 
photons in the process of photodissociation. Therefore, during 

computations of CO2 pj^otolysis we integrated up to X = 2050 S. 

If we consider the entire cross-section of attenuation due to 


photodissociation, then the rate of photolysis,at the surface of 

Mars is an order of a magnitude greater. 

In computing ^ (in an effective spectral interval 

UU 2 ^2 112^ 

X $ 2050 S), in the. sum, found in the index of the exponent ’.in formula 

(2), only the member whichetakes COzinto consideration is -essential (in view 

Qf its- predominance) .and' the- integral ’al-'opg Z; is‘ replaced by ~ °C02 

(Z) (Z)secX ( since CO2 is distributed according to the barometric 

UU2 

formula). Thus., it seems that these members are the functions of the 

number of CO2 molecules in the R column. These functions can be solved 

CU2 

beforehand, having been eliminated from integration for X and for Z 
in each step of the computation. On the other hand, when calculating 
Lq , Lg Q , Ljjq (in the spectral interval X < 2000 S, where the 
cross-sections of these components are much larger than the CO2 /8. 

cross-sections and are comparable among themselves), all of them are 
included in the sum in the index of the exponent; in addition, the 


distributions of these components differ from the barometric formula; 
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therefore, when computing these members one has to integrate both 
for Z and for X. 

To determine the average global distributions of 
we set the average daily profiles T(Z) according to [3/ 33] and 
the average values for the photodissociation rates L^, for which we. 

efr 

solve it when X=60° and we take \ the resulting value to find out that 

A 

the source is acting for one half the time of day, i.e. 

examination of the previous models (including 
models for the composition of the Earth’s upper atmosphere) and of our 
model of the thermal conditions of the Martian thermosphere shows that 
the average distributions found in this way are close to the conditions 
near the terminator in the middle latitudes. 

The members describing the binary chemi'cal reactions have the 
form: K^n^n^, and the triple reactions: - K^n^^n^n. The accepted 

system of the reactions and the coefficients of the reaction speeds, 
chosen from the latest laboratory measurements and evaluations [‘45-62], 
are presented in table 1, which demonstrates the majority of 
changes sharply depending on the temperature. 

The model includes 12 neutral components: CO2, O2, CO, O2 , 

0(^P), 0(^D), H, H2, H2 O, OH, HO2, H2O2, and also the ions CO2 , 

CO2H . And for 11 components the system of equations (l) was solved, 
and the concentration of CO2 was found from the barometric formula 
(Z = 0) and T(Z), according to the model [3, 33]. Evaluations 

showed that this is a very close approximation, because, in the 

« 

first place, the characteristic time for existence of CO2 in the 
Martian atmosphere is very large, many orders of magnitude larger 



than the characteristic time pro'vided. in the prohlem, and, in the 


/9. 

second place, the CO2 concentration is much larger than the concent ... 
trations of other components in almost all areas, and at high altitudes 
where the 0 concentration approximates that of CO2, the diffusion time 
is so small that n^^ is determined by its concentration in the lower 
layers, so that the barometric formula for it is not impaired for 
all practical purposes . 

The ion reactions [34, 35] (table l) significantly affect "the 
concentration of the neutral hydrogen n„ at altitudes greater than 

ll 

'^100 km. The CO2 concentration was taken close to the empiric profiles, 
according to the measurements of the "Viking" spacecraft [12] . 

The expressions for the diffusive streams in a multi-component 
mixhure of gases are very complex [63, 64] , however, the predominance 
of 0O2 in the Martian atmosphere allows us to use an approximation of 
the small component, hy using - the diffusion coefficients of the 
k- component in CO2 is very well approximated in the form 

= AT^n The expressions D^' for the components in the computations 
are taken from experimental data [65, 66] , and for several molecules 
for which there are no data, they are taken from the molecular 
characteristics 5 presented in '[63, 64.]; moreover, for these 
molecules of the short-lived components, D is completely insignificant 
(in view of the presence of photochemical balance = 1^)* 

The turbulent mixing coefficient K was given as an altitude 
function in the following manner: we considered that near the 

homopause the values of K can change from Ki=9-.10? to S.lO^cm^sec 
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according to the estimates [67] ; then the altitude of the hcmopause 
changes from Z 1=120-2^4-0 km; at the planet's surface the values for 
K can change from Ko=10^-10®cm^sec ^ [68-69]; from the surface to the 
homopause we considered that log K increases linearly (see figures 1-^ 
and table 2 ) . 

In computing the models we assigned the following boundary 

conditions: on the lower boundary (at the planet's surface, where 

Z = O) we assigned the absence of streams for all components, except /lO. 

H2,0. For water vapors we assigned a definite value of the concentration 

n„ ~(0), starting from the assigned humidity which, for different vari- 
1120 

at ions of the model, changes from 0.6 to TOU of precipitated water 
in accordance with the experimental data [16, 17] . At the upper 
boundary (Z = 250 km) we assigned the condition of diffusive equilibrium 
for all components, except H, Ha and 0. For H and Ha we used the 
departing stream because of the themal dissipation, according to 


the Jean formula: 



-.-totnjAL page 

^rPWR ftUALlTV 

Xc)eXf(-7c) f]H,H;, ^ ; (3.) 


where K is the Boltzman constant, subscript c is the value of the 

- —rj^- 

exobase ‘ ' /A ' ^ TT • ■ oxygen we used the departing stream. 


exobase ‘ ' /A ' ^ TT" • ■ oxygen we used the departing stream, 

equal to half the departing stream of the atomic hydrogen, i.e. 
:‘Po=|-(^k-+ 2 accordance with the McElroy conclusion 

on the presence of nonthermal escape of oxygen as a result of the 
fact that during the dissociation-recombination reactions 

® -^'CQ -4*. 0 H jV 0; 0., the O atoms obtain a 
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speed larger than the critical speed [30, 70] . The estimates of the 

number of atoms -with velocities higher than the critical, forming in the 

dissociation-re’combination reactions, ,are also close to the ^ value. 

o 

The choice of the initial conditions of n and n„ is made 

CO O2 

according to empirical data and is based on the following considerations. 

The dominant components of the Martian atmosphere enter into reactions 

as a result of which (ignoring several secondary reactions) can be written 

as 200^ 2C0 + 0^, . i.e. for each O 2 moleci^e, two CO . , ^ 

molecxiles are formed and disappear or 2AN^- and ' m 

. r'_. ‘ ■ c::: . . r . ■■ ..-i' CO ^ . O 2 

consequently, w . a- + A , where A is ahconstant.?) If / 

tU “ 'U2 

A = 0 (i.e. at the initial moment the Martian atmosphere would have consisted /ll. 
of pure CO^), then it would always be that ^ 2. However, the 

experiments show that in the Martian atmosphere the value of this rela- 
tionship lies within the limits of 0- 5-1-0. As initial conditions for CO 
and 0^ we also chose the latter profiles, so that constant A-= -2.12'*10^°cm ^ 

(a certain exicess of 0 atoms is assigned). In this case with = 

=2.12*10^‘’cm ^ we get = 2.12»12^®cm"^ and = 1, which is 

O 2 - 4* - - UO 2 U 

close to the experimental data. The computations showed that in all 
variations of the model the relationship + A) ^ 1 with’ good 

accuracy Cthe deviations do not exceed h%, see table 2), This shows 
that at first this difference in the inflow and outflow of water 
does not lead to a significant accumulation of 0? in the Martian ■ . 2 
atmosphere and, secondly - the computation scheme is conservative. 


tf) 


The relationship 2AC=A0 is a little more accurately fulfilled (locally) , 
which may be verified by writing their equality thusly: 


2f^ + 3f„ + f. 


0'^ 


0 - 


H 2 O 


+ 2f. 


H 2 O 2 


+ f . 


0 


^OH 


^^H02 - ^O(^D) 


- f 


CO 


= 0 


land substituting the expression for However, even this 

equality is impaired (integrally) because of the differences betw 


between 
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2. Ntunerical Method 


For a numerical solution the initial system of eg.uations 


(l) is •written as: 

■ - il. [f 4- ^Vi — / -ilK Vi’ i p V7jf "• , w 


where yjj.(2) is the steady-state solution of the boundary problem with 
zero streams on the border, with arbitrary initial conditions 
and zero and 

.1 o f ifaw. 1 )kH+khk1j4| 


The constant C is determined by selected initial, conditions. 

The expression (?) is a generalization of the barometric formula for 
the k-component. Writing the equation in the form (it) permits the 
isolation of the exponential dependence of the concentration on the 
altitude, which significantly improves the differential approximation 
of the initial system of differential equations, because many 
components of the Martian atmosphere are distributed according to the 
barometric formilLas O 2 , H 2 O, CO, H 2 at all intervals of the altitudes, 
K and 0 in the upper layer of the atmosphere). 

A stationary solution is sought by the determination method, i.e. 
we are integrating for t until we find a stationary solution. 
Differential equations are replaced by the differences on a net'work 


the amoxmt of 0 atoms leaving the atmosphere. Precisely because of this 
there is an excess of O 2 in the Martian atmosphere, described 
herein with the constant A. 
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of uneven altitudes (with thickening at the "bottom) . In most of the 
con5)utations we used the following network; = 0, 0.1, 0.3, 0.93 
2.T, 12.5, i8.T, 25, 32, 39, hj , 6o, 7^, 91.2, 105,. 120,- lUo., 

l60, 200, 250 Such' a concentration of points is made for 

improving the approximation of the system of diffusion equations 
by the differences at a set number of points and for decreasing theV 
possible error because of the poor approximation of the lower 
limiting border -condition (zero stream). 

We' use a^-conservatdve, unclear differential scheme of the second 
order' of .accuracy for Z: the, diffusion member is written along three 

points corre„sponding to [71] , but the nonlinear members are aligned in 
the middle point and at each intermittent step they are improved by 
iterations according to the Weirton method. The. solution of linear 
equations at'each iteration is conducted jointly by the exclusion method 
of 0auss. Thus, the differential scheme takes the form: 


p-fi, 






p 

P 




At . p 




rj3|K 
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here {Z ,m =’0.1,...M, Z_ = 0, = 250 km, AZm = Zm - Zm-l} is an 

m • 0 W 

*> 

■uneven network. 
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m is the nuinber of the point hy the altitude; A2m is the 
step by altitude; 

V is the number of steps in time; At is the step in time; 
. P is the number of iteration by nonlinearity; 

K is the number of the con^tonent; 


= OS' y 

Q r: - Q L ^Kjw»4, iSkJ 



1' 




^ CXma- 


m-i/£ 


I 


RisSi 

3/<J. 






m 
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— the expression for diffusive streams . 

To obtain a conservative differential scheme we must also write 
out the differential expressions for diffusive streams on the 
boundaries. The following expressions maintain the conservatism of the 
s cheme : 
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- on. the upper "boundary; 


0 

I'stj 




ji g 


on the lower "boundary. 






" S'M:vi=T^h(2c')-+'?h2.(^'^) _ 

The "boundary condition for H 2 O at the lower "boundary is written 

as: n„ _ ^ = Ci, where Ci is the assigned concentration. 

1I2U5 U 

In view of the large scatter of characteristic times 
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the linear differential system of equations, approximating the 
system of differential equations, is poorly stipulated. The 
poor stipulation leads to an increase in errors in rounding off with 
the growth of At. In conjunction with this at large times (larger 
than approximately 10^ sec) we had to compute with dou"ble accuracy 
(Fortran BESM-6, twenty- four decimal places). Since the num"ber of 
required operations for solving the system of equations according to 
the Gauss method is proportional to the cute of the number of 
equations, the multiple components of the problem also lead to 
serious difficulties. 

To overcome the difficulties we used the following scheme of 
computations : 

1. From the initial conditions we make one step in time AtJ^O.l sec.; 





if the iterations coincide, At is increased "by a factor of 1.5; 
if the iterations do not coincide (since at the given At sharp changes 


in . take place), the step is decreased hy 3 times; these 

computations are conducted! >on a coarse network - ten points along the 
altitude. 


2. With the appearance of large errors in rounding off 

(they appear "because at the given At the concentrations n^ change very 
little which is confirmed by the integrating at the point of time 
-At with a small step At; the iterations nevertheless do not coincide - 
the significant numbers n^ "dangle”) the computations are continued with 
double accuracy on the BESM-6 computer the speed of the computation is 
decreased by 4-5 times ) ; the computations are also continued on the 
coarse network. 

3. After obtaining the steady solution on the coarse network >: 
the number of points increases until we obtain the necessary accuracy 
(21 points were sufficient to arrive at accuracy). 

Receiving a steady solution with’ this -set ofhcomputations takes 
1 hour. 

After this we vary T, K - the time for solving different 

Jn.2^ 

variations is 6 minutes (time for first iteration at the 21st point 
on the altitude with double accuracy is '^1.5 minutes). 

In conclusion, we can say that the numerical method was 
practical in solving the problems with the number of components in 
the order of ten. This method is the strictest, does not use the 


/15. 


specifics of the chemical and diffusion processes and therefore 
is most general. Of course, other methods of solution are possible 
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•which use the specifics of the processes and therefore can he more eco’- 
nomical. For example, in the work [ 73 ‘, 74 ] the component cleavage 
of the chemical processes is used and the correction of the solution 
takes place hy examining the groups of active components, the concen- 
tration of which slowly changes in time. The work [ 72 ] uses the cleavage 
of the physical processes, i.e. the diffusion is separated from the ■ 
chemistry. After this, the diffusion eq.uations are cleaved hy 
component and are solved hy suhseq^uent scalar trial runs, and the 
equations of chemical kinetics are solved jointly. 

3 . Results and Discussion 

The altitude distribution of the concentrations n^ in 
several basic variations of the model are presented in figures 1-U 
and, in addition, the basic characteristic of several of the variants 
are presented in table 2, where the available experimental data are 
presented. The figxires show that the long-lived components {CO2, CO, /16. 
O2, H2, H2O) are distributed in the altitude according to the 
barometric formula, -whereas in the homosphere the altitude scale is the 
same for all of the components, and in the heterosphere - they are 
different; in the water vapor in the troposphere the altitude scale is 
determined by the curve of saturation and is approximately equal to 2 km, 
which is significantly less than H-. The remaining short-lived components 
have more complicated profiles, and the profiles 0 (^P), 0 (^D), H and OH 
have maximimis at the planet’s stirface. [C02].is comparable with 
[ 0 ] at 2=250 km. (in the basic variation). 
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2 . Let us compare the characteristics of our model with the 
experimental data which, as we showed in the introduction, are quite 
numerous and varied. As we see from tahle 2 , the values of the total 
content of O2 and CO in most variations of the model are within the 
limits of the distribution of the experimental data, with the 
exception of low-humidity variation (0.6‘y of precipitated water), 
which yields somewhat larger values. The values of the relative 
concentrations [0]/[C02] and [CO]/[C02l on the ionospheric maximum are 
also within the distribution field of the experimental data except'' 
for the variant with the minimal turbulent mixing coefficient in the 
thermosphere (K^ = 9 - lO^cm^/sec ) , •which yields values that are several 
times larger. It is entirely possible that there variants may be reap, 
lized also at some time. 

The values for the concentration of atomic hydrogen on the upper 
boundary at an increased temperature and decreased humidity 
approximate the experimental values = (2 - 3)*10‘*cra but in most 

of the variants they exceed these values several times. The basic 
model satisfies the experimental data best of all (the average for 
K in the upper atmosphere and the minimal in the lower atmosphere). 

In this model all of the concentrations lie within the limits of /IT * 

the scatter zone of the experimental data, except a slightly increased 
content of O3 (T.^ y-atm. ) and a slightly increased “ 3 - 7 * 10 “* cm 

The ozone concentrations in our model change from 21 y-atm. at 
low humility (equal to 0.6 y of precipitated water) to 4.7 y-atm. 

O3 at 15 y or more of precipitated water (figure 5 )- These values 
are within the scatter limits of the experimental data ( 3-60 y-atm. O3). 
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True, the Talues at high humidity (approximately 5 p-atm.) are 
slightly higher than the experimental (<3 U-atiri) , however, it 
was found that complete agreement is obtained if we take a warmer 
troposphere; this is completely cut, i,e. in our model we took 'the 
global, average value at' the surface of the planet T^ = 220 K, 
whereas at the lower latitudes where there is a lot of humidity 
[16, 17], significantly higher temperatures = 250 - 2T0°K are usual 
[3, 5, 75] . The ozone content at the increased temperature increases 
because of the reverse temperature dependence of the speed of the 

reaction K... . 

11 

This permits us to conclude that the data of our model are in good 
agreement with all oT the available experimental data. 

3. The model shows that the variations in intensity of the 
turbulent mixing in the thermosphere within the possible limits 
(K^ = 9.10^-8.10^cm^sec result in noticeable variations of the 
concentrations of several components (compare the variations 2 and 3 
in table ,2:).. The 0 and CO concentrztions decrease especially strongly 
with an increase ih' K (approximately by 6 and 3 times, respectively). 
Apparently this takes place as a result of an increase in the transfer 
of the CO 2 molecules upwards where they are subjected to 
photodissociation, and a transfer of CO and 0 downward where they 
associate as a result of the cycles of chemical reactions described 
below. Lower values of the H concentration at the exobase also 
correspond to a greater intensity of the turbulent mixing in the 
thermosphere . 

The turbulence in the lower atmosphere also plays an important 


/18. 


18. 



role.*. With the reduction of K the total content of CO, O2, and O3 

increases and the level of H2 declines, ^whereas the H and H2 at 

the exohase decrease. In our model vhen K = 10 cm sec •, in the 

o 

troposphere the best concurrence "with the experimental data is obtained. 

If. Studies of the effect of the variations of temperature in the 
upper thermosphere (the average global values of "which can vary 
from T =262°ICat low solar activity to =35,9°Kat high solar activity 
[76]) shows (compare variants U and 5 of the model) that with an 
increase in within the prescribed limitations, the concentrations 
of several components at high altitudes increase (for example, 

[COal^^Q increases by approximately 5 times) because of an increase 

in the altitude scale, and the concentrations of the atomic and molecular 

hydrogen noticeably decrease (by 5 and 1.6 times, respectively) 

because of an increase in the stream of the thermal dissipation, the effect 

of which prevails over the effect of ah increase in the altitude 

s cale . 

5. The effect of the increase in the humidity is manifested, 

as we noted above, by a strong decrease in the ozone content'* . ' 

(figure 5). In addition, "with an increase in the hvutiidity from 0.6 p 
to 35 y of the precipitation, the level of H, H2 and OH significantly 
increases; thus, the concentrations at the upper boundary 

[H2]2 ^q increase approximately T times, and the concentrations at 
the lower boundary [OH] ^ increase 10 times. 

6. We studied the effect of the variations of the ion compo- 
sit ion, or the amount of [CO2] , on the concentration of the atomic 
hydrogen in the thermosphere. Since the hydrogen in the thermosphere 
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is basically formed in the reactions (3^, 35 )s then with a decrease 

in [cot] ty 3 times, the concentration of decreases hy 2 times /19« 

("variant s 9 and 10 ) . 

7. The model. shows that the effect of the reaction of the 
photodissociation predominates at altitudes higher than approximately 
50 km, hut in velocities the chemical relations are comparable with 
them at lower altitudes (figure 6). The difference between the 
speeds of fomation and destruction of the components is compensated 
for by the diffusion streams below 120- 140 km - as a result of the '■ 
turbulent diffusion), which, as we see from figure 7, transfer the 
CO and 0 downward; O2 is transferred upward at altitudes more than 
'^'60 km and do"wnward at' lower, altitudes ; H - upward at altitudes more 
than '^105 km and downward at lower altitudes ; H2 - upward at altitudes 
more than approximately 20 km and downward at lower altitudes; O3 - 
upward at altitudes more than 40 km and downward at lower altitudes. 

Some of the reactions introduced hy us turned out to he insig- 
nificant (Ts 8, 12, 16, 21, .27, 29, 33)5 meaning that their contribution 
to the balance of any component in the entire range of altitudes 
examined is less.ithan 1^. These reactions may he excluded from 
further computations . " ‘ v 

In previous reports [2, 29 - 32] , we indicated several catalytic 
cycles aiding the association of CO and 0, specifically: 

1 st cycle-: 

c H + O 2 + M HO 2 + M 
HO 2 + 0 -s- OH + O 2 
CO OH ^ CO 2 + H 
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resulting: CO + 0 CO 2 

then the II cycle: 

H + O 2 + M ^ HO 2 + M 

HO2 + HO2 H2O2 + O2 

H 2 O 2 + hV 20H 

CO + OH CO 2 + H 


• . ('.resulting: 2C0 = 02 ;^ 2 CO 2 

and, finally, the III cycle: 

0 + O 2 + M — O 3 + M 
O 3 + H - OH + O 2 
CO + OH - CO 2 + H 


resulting: CO + 0 - CO 2 

This yielded opposing opinions about the relative humidity of 
these cycles. Our model shows that all three cycles play a definite 
roles, and at the surface in the 2nd cycle 63^ of the OH molecules 
form, in the 1st cycle - 35^ and in the 3rd - approximately 1^; at 
an altitude of several hm the prevailing contribution is given by the 
1st cycle, the 2nd cycle becomes negligible; at an altitude of 
approxiamtely 30 km the 3rd cycle '-s contribution is comparable 
to that of the 1st cycle and at higher altitudes (up to 100 Ion) the 
3rd cycle prevails. 

An important role in the chemistry of the Martian atmosphere 
is played by excited oxygen 0(^D) atoms in spite of the fact that 
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their concentration is very small (ahoat 1 cm ^ ) . These atoms are 
formed in the photodissociation of O3 (reaction 3) vith a speed 
approximately of l.l’lO^cm ^sec and immediately, approximately 
with the same speed, they are destroyed in the impact deactivation with 
CO2 (No. 18), however, approximately S'lO"* atoms cm ^sec enetring 
reaction 19 effectiveljr 'break up the hydrogen molecules and approximately 
3.5’10^ atoms cm ^sec~^ as a result of reaction 20 break up the H2'’ 


molecule . 

8. The global, daily-average streams of dissociating radiation 

which we used can increase the concentrations of the short-lived 

components. Therefore, we computed a variation (figure 4-, No. 11 in 

table 2 ) , in which we used the average daily , L ^ and L 

^3 , -^ 2^2 ii'Jz 

(i.e,”L^ (x = 60 °)), twice greater than the daily average. In this 

variation we obtained an ozone layer (N^ =3.^ p-atm. ) close to the 

O3 

experimental level. It is also interesting that in this variation 
(in comparison with variation No. l) that the concentration of CO at 
the surface decreased two times, and the concentration of O2 decreased 


1.5 times, whereas in the thermosphere the concentration of all of these 


components practically does not change. Apparently, this decrease 

took place as a result of the cycles of chemical reactions mentioned 

above. Further improvement of this factor is possible onyl within 

the framework of a model describing the daily variations. 

To compare our model with that of Liu and Donahue [32] we, computed variatioi 

(no. 12 in table 2) with parameters close to those used in [32] 

(K = U*10®cm^sec \ N„ - = 15P‘ of the precipitated water, close profiles 

n2^ 

of T(Z), (z)). The concentration of most of the components 

OU2 

in this variation satisfactorily agree with those in [32] , and the 


22 . 



ratios of ll|.0^ lltO almost one order less 

than that of [32]. The concentration of [CO 2 ] = 1.65*10 ^['COz]^ 
and IO 2 ] = 6.10”** [COzl^ are considerably lower than in [32] apparently 
because of lesser speed of the photedissociation in our model, 
equal to l.T*10^^cm ^sec~^ (see the speeds of the reactions in figure 5). 
In addition, it can be seen that in our model some additional 
conditions were not applied, which were used in [32] . 

10. Measurement s using the mass-spectrometer on the Viking 
spacecraft showed that present ..in. the "^Martian atmosphere 'are nitrogen 

•'^[^ 2 ] ^ 2-6^ [CO 2 ]) and nitrous oxide ([WO] lO^cra ^ at an altitude 
of 110 kilometers ) . Evaluations showed that these small concentrations 
of WO do not significantly alter the picture obtained from our model. 
However, the concentrations of several small components (particularly 
“ HO 2 ) can change somewhat. Moreover, there is a possible effect of 

“f* 

WO on the ionosphere through the reactions W 2 + 0 ->■ WO + W and 

W 2 + HO ^ ho'*' + HO 2 which have high speeds in the order of 10 ^°cm ^sec ’ 

11. For different variations of the described model we 
performed evaluations of the departing stream of oxygen atoms 
arising with the dissociative recombination of the O 2 and CO 2 ions 
(an exact computation is possible only with an exact determination of 
the distribution of the ion concentration and the electron and 

ion temperatures). It turned out that the is much less variable 
than the departing stream of hydrogen atoms ($g+2§jj^) (because the 
first one is determined ultimately by the concentration of the basic 
component CO 2 and by the speed of the ion formation, whereas the 
second one depends more' on the neutral composition and its variations). 
In several variations the ratio of these streams is entirely 
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fo/($rr+2$._ ) 0.2 — 0.3, that is less than 0.3, "which was permitted 

n il 2 ^ 

earlier. It is possible hecanse of this consideration over a long 
period of time to build up an excess of oxygen atoms, noted above (see 
page ll).* 

In conclusion we would like to express our thanks to A. I. 
Poroykova (institute of Chemical Physics, Academy. of. Sciences USSE) , 
for help in selecting the coefficients of the speeds of the chemical 

I 

reactions, E.E. Shnol (institute of Applied Mathematics, Academy 
of sciences USSR) for his consultation on the numerical methods and 
T.M. Donahue (University of Michigan) for familiarizing us with h'is 
report [32] before its publications. 
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Table 1 


■ 

• 




No. 

Reaction 

Speed Reference (•*) 

I. 

C02+lii)C^O~2062) CO-j O 

X,» 3,62.10-'^ [37;38„39l. 


O 2 +W (^0-^2062) OtO 

2 , 00 , 10 "® [36,40,4X3 

3/ 

O 3 +KI (2000-3600) ->-■ 0(^D )+02 

4,12.10"^ 

[36] 

4. 

H20+h9 (1515-2000) -s. OH+H 

2,06.10"® 

[ 42 I 

5. 

H202+f)9(2000-3600) OH+OH 

7,41.10"^ 

P'3] 


H02+W( 2000-3600) O+OH 

5,65.10"^ ■ 


7. 

CO 'f 0 + OO 2 ^ CO 2 + ^^2 

1 , 01 . 10 " 

'[45,461. 

8 „. 

CO + 0 CO 2 + 

6,64.I0"^®eX/'<"^f^) 

&7] 

9, 

CO -}• OH CO 2 H 

(-^®) 

M 

10 . 

0 + 0 i- C 02 “^ O 2 + CO 2 

3.10“*^^ ( j 

[ 49 . 50 ] 

II. 

0 ©2 "J" CO 2 -fr- O 3 + ^^2 

0 . 300t . 

1,4. ( ‘ 

300 

5[49,5qj 

12 . 

O'}- 0 3 2 O 2 


[5ll 

13. 

0 ^ H 2 -^ OH + H 

7,0.I0-Il£X/^(-5000) 

[30] 


0 4 OH O 2 + H 

4,0.10"*^^ 

[ 52 . 53] 

15. 

0 ^ HO 2 OH ■{' 0 -^ 

8,0.I0"II^)f^(-522) 

[3] 

16. 

0 + ^ 2^2 ^ 2 ^ ^2 

’ 

2,3.I0-n^(-32J0) 

l55] 

I?; 

. 

0 + H 2 O 2 OH + HO 2 • 

p <; ■* 6 32CK)% 

?,3.I0 -vgjj^C,—) 

[ 35 ] 

18. 

0(%) ■)+ C02-i-0(%)+ COg 


[36] : 

19. 

0(^D ) +E 2 O 20H 

!3,0.I0-^° 

j 

[ 57 , 5 ^ 

' 20 . 

0(^ ) + H 

^3,0.10”^° 

[57,58] 
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Table 1 'continued: 



*) For J - a reference to the source from which are taken to compute 
' T - the value of X at the boundary of the Martian atmosphere. 

• bo' S: 
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I' Variation T3o, . 


T , “K 




WhoO^ p -5 

JV 02 Ccn“2). 10-20 

iVo^/ZVco.-m*^ 

W 00^02*®^^ 

MoS.v -a™) 




2iI40’ 


[h]25oCch“^).IC-^ 

[H2“]25o(cm'*^,-IO-5 

[ 03 loCcr^),IO~^ 

:[H2023oCcr^)oIO*’5 

if 03 man ■(CK*'^)«I0“'^^ 

.[0Hlm^5cu-“^„IO-'5 

;[0(% (CK^) 

.Wco'^llo^ ' 'T 
'NcoA^ Mq^ -^ ^') •- 


15 

3,13 


I 

15 

2.43 

1.2 


2,74 


1.3 

4.7 

1.3 

1.1 
7,9 

5.2 

3.7 

5.0 

3.8 

6.1 
4,7 
2.0 
2,1 


322 

0,9 

I 

15 

2,52 

1,2 

2,91 

Io4 

4-,7 

6o3 

3.2 

9.1 
5.5 

3.4 

5.1 
3,9 

6.5 
7,5. 
1,7 
2„2 


I 

15 

2,43 


I,X I 1.2- 
IoQ’0!r499-7 


5.0 
3.8 

6.1 
4.7 
2,0 
2,1 

I, I 

ijd 










Talile 2 












Figure 3. The same as in Figure 1, with minimal moisiiure- - 
of HgO precipitation (variation No. 4 in table 2). 
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■ y 

Figure 5 . Dependence of the tothl,.* ozone content on the moisinre. 

The solid curve corresponds to the coefficient of the turbulent 
mixing in the troposphere K =10^cm^sec~S + ~ K^=10^cm^sec~^ ; 

X is the average daily speeS of photodissociation of O3, h^Oo 
and HO2 at KQ=10^cm^sec“^ ^ 
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Figure 


Log V, cm^sec 


6. The speed of the basic reactions (variation No. 2 in table^2) 
Fotatibns: + - the speed of the first reaction. 
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